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ABSTRACT 


The Helmholtz theorem is rederived with rigorous vector analysis The theorem is 
valid everywhere within any arbitrary mathematical boundary. Applications of the 
theorem to hydrodynamics are discussed. 
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INTRODUCTION 


One of the theorems of Helmholtz states that a vector field with moderate restrictions 
on continuity can be expressed as the sum of the gradient of a scalar potential and 
the curl of a vector potential. This theorem has been known for a long time, and is 
derived in texts 1-3 on vector analysis An elegant proof was prepared as part of a 
course which was taught in 1958 at the Naval Proving Ground Publication of the proof 
at this time is justified by the impending application of the Helmholtz theorem to the 
computation of flow around a ship. 

The existence of the scalar and vector potentials is proven by the derivation of 
formulae which express the potentials. Only vector identities are used in the derivation. 
The theorem is valid within any closed mathematical boundary. 

The algebra and the calculus of scalars, vectors, and tensors are the subjects of 
many texts Proofs are both geometrical constructions and component manipulations. 
The most relevant theorems which are background for the derivation are presented 
herewith Formulations are stated in Gibbs notation. Symbolic expressions are invariant 
with respect to the choice of coordinate system. 

RIEMANN INTEGRATION 

Let a number of variables vary over a range of integration. Let the range be divided 
into elements such that the increment of each variable in any element is less than c. 
An analytic function of the variables can be approximated within each element by a 
Taylor series expansion. The true integral of the function within each element differs 
from the trapezoidal integral by infinitesimals of higher order than e. The true integral 
reduces to the trapezoidal integral in the limit as c goes to zero. The structure of 
each element is immaterial as long as the elements fill the range of integration. This 
is known as Riemann integration 4 . 

VECTORS 

A vector is a quantity with magnitude and direction. The magnitude and the direction 
are invariants of space. 

The scalar product of two vectors a and b is given by the equation 

a b = |a!|b| cos 9 (1) 

where |a| and |b| are the magnitudes of the vectors and 9 is the angle between them. 

The vector product of two vectors a and b is given by the equation 

axb = n |a|'b| sin 9 (2) 

where n is a unit vector normal to the plane of a and b, |a| and |b| are the magnitudes 
of the vectors, and 9 is the angle between them. The direction of n is such that the 
direction of b is obtained from the direction of a by a right-handed rotation about 
the vector n. 

The scalar-vector product of three vectors is given by the equation 

a bxc = (a n) |b||c| sin 0 (3) 

where n is a unit vector normal to the plane of b and c, while 9 is the angle between 
b and c. The scalar product a n is the separation of two parallelograms of area bxc. 
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Thus the scalar-vector product abxc is the volume of a parallelepiped with edges a. b, c. 
The triple vector product of three vectors is given by the equation 

ax(bxc) = (axn) b||c| sin 6 (4) 

where n is a unit vector normal to the plane of b and c. while 6 is the angle between 

b and c. The vector product axn has a magnitude equal to the projection of a on the 

plane of b and c, while it has a direction perpendicular to the projection of a on the 

plane of b and c. Planes which contain the tip of a and are perpendicular to b and c 

intersect b and c at the projections of a on the directions of b and c. Lines which 
contain the origin of a and are parallel to the planes complete a parallelogram whose 
sides are the components of a in directions perpendicular to b and c After rotation 
through a right angle the components of a become the components of a> n in directions 
parallel to b and c. The triple vector product is given by the equation 

ax(bxc) = ba c - ca b (5) 

which resolves the triple vector product into components in the plane of b and c. 

A tensor is an operator which converts one set of vectors into another set in 
accordance with a linear transformation. The tensor has a matrix whose elements are 
the coefficients in the linear transformation. Typical tensors are rotations and 
deformations. 

The tensor product of three vectors is given by the equation 

ab c = (c b) a (6) 

where c is transformed into a vector in the direction of a by a special tensor with a 
matrix whose elements are the products of the components of a and b. Any letter 
could be used to represent the tensor, but the dyadic notation ab is more useful. 

Let i, j, k be orthogonal unit vectors in the directions of increasing coordinates 
x,y, z in a right-handed Cartesian coordinate system. Then the vectors have the 
properties 

ii = jj = kk=l ij = jk = ki = 0 (7) 

ixj = jxj = k*k = 0 ixj = k jxk = i kxi = j (8) 

Any vector can be resolved into components along the unit vectors. 

All nine pairs of vectors among i. j.k form a fundamental set of unit dyadics. Any 
tensor can be resolved into components along the unit dyadics. The identity tensor I 
is defined by the equation 

I = ii + jj + kk (9) 

The identity tensor transforms every vector into itself. 


GRADIENTS 


The gradient is an operator which gives the differential change of a function for a 
differential displacement in space. The gradient of a function A is expressed by the 
equation 


dk .3A , 9A 
VA = i—+ j— + k — 
dx dy dz 

where the function A may be a scalar, vector, or tensor. 


(10) 
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The divergence of a function A is expressed by the equation 


. dA . 3A , dA 
VA = i — + j — + k — 
dx dy 02 

where the function A may be a vector or tensor. 

The curl of a function A is expressed by the equation 


VxA 


OA . 3A . dA 
i x t~ + J*— + kx — 
dx dy dz 


( 12 ) 


where the function A may be a vector or tensor. 

The partial derivatives of a continuous function are independent of the order of 
differentiation. Thus the function A satisfies the identities VxVA = 0 and V-VxA = 0. 

The Laplacian of a function A is expressed by the equation 


V-VA 


d*A d z A a 2 A 
dx z dy z dz* 


(13) 


where the function A may be a scalar, vector, or tensor. 

Laplace’s equation is V-VA = 0. An important solution of Laplace’s equation is given 
by the equation 


1_ 1 

r + y 2 + 2 Z 

which may be verified by substitution in the Laplacian. 


(14) 


GAUSS THEOREM 


Let a volume of space be enclosed by a boundary surface. Integration of the 
divergence throughout the volume is expressed by the equation 

J'A* T = ///.£ ***»* + J J J J-fJ dxdy ’“ + J J j V fz ***** <15) 

where dr is a differential element of volume. Partial integration leads to the equation 

Jv-Adr = J J* [i'A j* > dyd 2 + J J [j-aJ > dzdi + J J [k-Aj > dzdy (16) 

where subscript (1) and superscript (2) indicate the difference between values of an 
integral at the opposite ends of the line of integration. The surface elements 


dydz 


dzdx 


dxdy 


(17) 


can be replaced in the surface integration by the parallelograms which are the 
projections of a surface element da on each coordinate plane. Thus 


/*•**--/ 


da A 


(18) 


where dr is the differential volume element and ds is the differential surface element. 
The vector ds is directed outward from the volume of integration. This is known as 
the Gauss theorem. 
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Let a surface be bounded by a closed circuit. Let u, v be right-handed coordinates 
in the surface. A vector surface element is defined by the equation 


ds 


\du dv} 


dudv 


(19) 


where r is the position vector of the element with coordinates u. v. Application of the 
triple vector product transformation leads to the equation 


/9r dr\ dr dr . . 9r 9A dr dA d / dr \ d (dr \ 

\9u dv) du dv dv du du dv du\dv J dv\du } 

Then integration is given by the equation 


f l 

* f dr 1 <2) / 

’ r ar 1 <2) , 

r dr f 

I ds-VxA = I 



b-A du + (t) 

Lav J (1) J 

lau J (1) j 

I du J 


( 20 ) 


( 21 ) 


where subscript (1) and superscript (2) indicate the difference between values of an 
integral at opposite ends of the line of integration. Thus 


J* ds-VxA = (j) dr-i 


( 22 ) 


where ds is the differential surface element and dr is the differential circuit 
displacement. Circuit integration is right-handed relative to the surface normal. This 
is known as the Stokes theorem. 


HELMHOLTZ THEOREM 


Let a volume be enclosed by a surface and let a point within the volume be enclosed 
by a sphere of radius e. In the space between the two surfaces a vector a satisfies 
the Green’s theorem which is expressed by the equation 


J V-(VU a) dT = J V-W ad t + J VU-Vm dr = J 


da-VU a 


(23) 


Let U be given by the equation 



V-W = 0 (24) 


where r is the distance of the point of integration from the center of the sphere. 
Within the sphere the integral 


* 



(r= e)(25) 


is the solid angle of the sphere. Therefore the value of a at the center of the sphere 
is given by the equation 


a 


_ 1 _ 

471 



•Va dr - 



(26) 
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Transformations of the integrals are made with the aid of the following identities. 

v|av(i)J = (Va)v(^) + aw(i) (27) 

v.j 1 „.v(i)j=v..v(i) + w(i).. (28) 

Inasmuch as W(£) is a symmetric tensor, it may be cancelled to give the equation 

J (Va)v|-j dr - f Va-v|-j dr = J dsaV^-j- j - J* ds a-V^ j (29) 


Transformations of the integrals are made with the aid of the following identities. 


7^—^x(axds) = aV^i j-ds - ds V^- j- 


(30) 

(31) 


|x(Vxa) = Va-V^j - V^J-Ca 
Therefore, the value of a at the center of the sphere is given by the equation 

8 - £ J v (;) V a dr - £ J v ( rh*"*■ - s J - -h J v(i)x<«xd») 02) 

Let the potentials <p and A be defined by the equations 


-in 

■in 


V-a dr 

Vxadr + 


1 f 1 . 

- — - a ds 

4n J r 

Lfi 

in J r 


axds 


(33) 

(34) 


where r is now the distance of the center of the sphere from the point of integration. 
Then the vector a is given by the equation 

a = - V<p + VxA (35) 

This is known as the Helmholtz theorem. 
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HYDRODYNAMICS 


Let a solid body be moving at constant speed U through a fluid which was initially 
at rest. A Cartesian coordinate system can be set up with origin fixed in the body. 
Let i be in the direction of motion, let j be in the direction to the right, and let k 
be in the direction downward. Then relative to this coordinate system the flow is 
backward. The flow consists of two parts. The first part is a uniform flow with a velocity 
which is constant everywhere, and the second part is a local flow with a velocity which 
diminishes rapidly with distance. Let there be a nonslip boundary condition at the 
surface of the body. Then the uniform flow has a velocity - Ui everywhere, while the 
local flow has a velocity + Ui at the surface of the body. The velocity v in the local 
flow satisfies the Helmholtz theorem 


v 


_1_ 

4n 


Mh dr 


_i_ 

4n 






(vxds) 


(36) 
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which gives the velocity at a field point in terms of velocities at all integration points. 
Volume integrals are evaluated throughout the fluid. Surface integrals are evaluated 
over the body, but vanish over the surface at infinity because of the diminution of 
velocity with distance. 

Insofar as the fluid is incompressible the divergence Vv is zero and the first integral 
vanishes everywhere. Next to the body there is a boundary layer. Within the boundary 
layer vorticity is diffusing outward. Outside the boundary layer the curl Vxv is zero 
and the second integral is limited to the boundary layer. 

In the limiting case of a thin boundary layer, let the volume element be a prism 
which straddles the boundary layer. Let the dimensions of the prism be e perpendicular 
to the boundary, d in the direction across the flow, and A in the direction along the 
flow. Then the volume element and the surface element are given by the equations 

dr = £<5A ds = <5An (37) 

where n is a unit vector in the inward direction. Application of the Stokes theorem 
to a rectangle of width e and length X leads to the equivalence 

Vxv dr = - vxds (38) 

For a field point far outside the boundary layer the second integral cancels the fourth 
integral, but for a field point inside the boundary layer the sign of V(i) is reversed 
and both integrals contribute to the velocity. The predominant integral outside the 
boundary layer is the third integral, which expresses the velocity as the integral over 
the body of radial flow from a continuous source distribution with the source density 

~ i n (39) 


where n is a unit vector directed outward from the body. Evaluation of the integrals 
is a quadrature which is complicated by the presence of an inverse square in the 
integrands. The inverse square introduces a spike in each integrand and the spike 
cannot be integrated by summation over finite intervals. A possible technique for 
integration through a spike will be the subject of a future report. 

The boundary layer flows into the wake behind the body. The structure of the wake 
depends upon the separation of flow in the boundary layer. The body may have a wave 
train. If the body is self-propelled, there is flow into the propeller. The uniform flow 
is disturbed by these local patterns of flow. The Helmholtz theorem may by applied 
to any or all of the patterns which disturb the uniform flow. 

A Rankine ovoid is that streamline of zero stream function which is generated by 
the ideal flow around a point source and a point sink. The NSRDC 5 made a model with 
the same profile as the Rankine ovoid and towed the model in the towing tank. The 
real model had a boundary layer and wake which would be absent from the ideal flow 
around the source and sink. 

The velocity v for the ideal flow is given by the equation 

v = - U\ - (40) 


where <p is the potential of the source and sink. The surface normal at the boundary . 
is given by the equation 


bxv 

|bxvj 


(41) 
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where b is the binormal transverse to the streamline. For the nonslip boundary 
condition with a thin boundary layer the real flux is just equal to the ideal flux as 
indicated by the equation 

l/i-n = - Vip-n (42) 

Thus the normal component of velocity is the same for both ideal flow and real flow. 
It is well known that two solutions of Laplace's equation which have the same normal 
gradients at a boundary surface are identical to within an additive constant. The real 
flow outside the boundary is the ideal flow from the source and sink. 

The NSRDC 5 measured the elevation in the wave train in the towing tank, and a 
computing program at NSWC 8 was used to compute the wave elevations for the point 
sources. The model was towed with a twisted cable and the boundary layer was fully 
turbulent. The forward motion in the boundary layer met partially the boundary 
conditions, and less source distribution was required to complete the boundary 
conditions. The experimental elevations were a few percent less than the computed 
elevations. Turbulence in the wake was revealed by photographs of smoke traces 7 . 

DISCUSSION 

Among the texts which have been inspected only two had derivations of the Helmholtz 
theorem, and they derived formulae for an infinite space, or for a field which was 
finite within a boundary and zero outside the boundary. In this report enough 
background material has been included to define notation and to clarify concepts. 

The gradient of velocity satisfies the identity 

Vv = |(Vv - V’v) + |(Vv + V’v) (43) 

where V’v is the transpose of Vv. The antisymmetric part is the rate of rotation and 
the symmetric part is the rate of strain. The differential change of velocity dv for a 
differential displacement dr is given by the equation 

dv = dr-Vv (44) 

Application of the triple vector product rule to the antisymmetric part gives the 
equation 

dr-|(Vv - V'v) = |(Vxv)xdr (45) 

Thus the angular velocity <•> is given by the equation 

u = |Vxv (46) 

It has become fashionable among hydrodynamicists to use omega for the curl itself, 
but this is a deplorable deviation from standard notation. The use of omega for angular 
velocity is pre-empted by its universal use in classical and quantum physics. It would 
be better to use gamma as defined by the equation 

y = Vxv (47) 

where y stands for the density of circulation. 

CONCLUSION 

The Helmholtz theorem is a useful law to which the flow in any flow field must 
conform within any mathematical boundary. 
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